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ABSTRACT 

Public data from the 2dF quasar survey (2QZ) and 2dF/SDSS LRG & QSO (2SLAQ), with their vast 
reservoirs of spectroscopically located and identified sources, afford us the chance to more accurately 
study their real space correlations in the hopes of identifying the physical processes that trigger 
quasar activity. We have used these two public databases to measure the projected cross correlation, 
uj p , between quasars and luminous red galaxies. We find the projected two-point correlation to have 
a fitted clustering radius of r , = 5.3 ± 0.6 and a slope, 7 = 1.83 ± 0.42 on scales from 0.7-27/i _1 Mpc. 

We attempt to understand this strong correlation by separating the LRG sample into 2 populations 
of blue and red galaxies. We measure at the cross correlation with each population. We find that these 
quasars have a stronger correlation amplitude with the bluer, more recently starforming population 
in our sample than the redder passively evolving population, which has a correlation that is much 
more noisy and seems to flatten on scales < 5ft _1 Mpc. We compare this result to published work 
on hierarchical models. The stronger correlation of bright quasars with LRGs that have undergone a 
recent burst of starformation suggests that the physical mechanisms that produce both activities are 
related and that minor mergers or tidal effects may be important triggers of bright quasar activity 
and/or that bright quasars are less highly biased than faint quasars. 
Subject headings: galaxies: general - large-scale structure of the universe - quasars: general 



1. INTRODUCTION 

The connection between quasars and their host galaxy 
evolution is notably demonstrated by the tight correla- 
tion of black hole (BH) mass to host bulge mass, the 
y^BH-P relation (e.g.,Gebhardt et al. 2000, Ferrarese & 
Mcrritt 2000). A question that emerges from this ob- 
servation is 'how does the BH at the center of a galaxy 
know what mass bulge it is sitting in?' A reasonable 
answer is that both of these masses grow together, and 
that the growth of both is triggered by the same mecha- 
nism. Since AGN are believed to be the sites where BHs 
are actively acquiring mass, we are motivated to look 
to these phenomena for clues to the physical mechanism 
that drives this correlation. 

A way to explain this connection of BH and bulge mass 
is to invoke major mergers of galaxies. Models of hier- 
archical growth of structures and galaxy evolution have 
suggested that quasars are fueled by galaxy major merg- 
ers (e.g. Haehnelt & Kauffmann 2000). During these 
events, gas is driven into the center of the galaxy, feeding 
the massive black hole within. Kauffmann & Haehnelt 
(2002) use such models to show that the quasar lumi- 
nosity function can qualitatively be reproduced, albeit 
with a very simple prescription for fueling and time- 
scales. In addition, recent models use the consequence 
of quasar/ AGN phenomena during galaxy mergers to al- 
leviate the long standing problem of an over-production 
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of high luminosity, blue galaxies while under-predicting 
the numbers of Milky Way-like galaxies in galaxy evolu- 
tion simulations (e.g.,Croton et al. 2006; Di Matteo et 
al. 2005). By providing a heating source for cool gas, 
quasars/ AGN added to these models are able to halt 
the build up of starformation in their hosts, arresting 
the growth of these luminous blue galaxies. The sugges- 
tion of a further connection between quasar/ AGN activ- 
ity and galaxy evolution is encouraging and leads to the 
expectation that quasars and AGN are found in over- 
dense regions where the merger fraction is expected to 
be high. However, there is also evidence that mergers 
operate more efficiently in overdense regions where the 
galaxy densities are not too high and so the encounter 
velocities are lower, like in galaxy groups or at the edges 
of rich clusters (e.g. Gnedin 2003). These are regions 
less likely to be occupied by red elliptical galaxies, more 
numerous toward the centers of rich clusters. Although 
in the most dense regions of clusters merging fractions 
are expected to be lower, the cumulative effect of weaker 
interactions like minor mergers or harassment may be 
important triggers of active nuclei. 

Recent large surveys like Sloan Digital Sky Sur- 
vey (SDSS), 2dF, 2QZ and 2df/SDSS LRG and QSO 
(2SLAQ) provide the opportunity to obtain good statis- 
tics on the environments of luminous and dim AGN. 
However results of some recent studies on the close envi- 
ronments of AGN provide a mixed picture. Kauffmann et 
al. (2003) use an SDSS sample to show that the hosts of 
narrow line [z < 0.3) are predominately massive, early- 
type galaxies (which are generally found in rich environ- 
ments). However, for the most luminous of these AGN, 
there are a high fraction in low density environments 
(Kauffmann et al. 2004; Coldwell et al. 2002). Further- 
more, Miller et al. (2003) find that the fraction of galax- 
ies that host AGN (z < 0.1, M r = -20) is independent 
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Fig. 1. — Histogram of redshifts for the 420 2QZ quasars (solid) 
and 4985 2SLAQ galaxies (dash-dot) used in this study. The me- 
dian redshifts for each sample are 0.59 and 0.55, respectively. 

of environment. While, Li et al. (2006) find that narrow- 
line AGN are more weakly clustered than a matched sam- 
ple of SDSS galaxies (0.01 < z < 0.3, -23 < M r < -17). 

On the more luminous end of the active nuclei range, 
Serber (2006) finds that luminous quasars (z < 0.4, 
Mi < — 22) cluster like L* galaxies for scales < IMpc. At 
high z, Coil et al. (2007) report that the SDSS quasars 
(0.7 < z < 1.4, -25 < M B < -22) that overlap with 
the DEEP2 survey are found to cluster more like blue 
galaxies than red. 

Here we explore one simple test of the major merger 
models of quasar triggering; Does the spatial distribution 
of quasars and LRGs support the picture where quasars 
eventually become LRGs? Our method is to measure 
the real spatial cross correlation of 2QZ quasars from 
the 2QZ survey with galaxies targeted by the 2SLAQ. 
The two main advantages of these large surveys in mea- 
suring correlation functions are 1) the large number of 
quasars and galaxies available in the datasets and 2) the 
availability of spectroscopic redshifts for a large fraction 
(> 90%) of the sources in the survey volume. Other stud- 
ies (e.g. Brown et al. 2000, Padmanabhan et al. 2008) 
rely on statistical photometric redshifts to determine the 
association of galaxies with quasars. 

In the next section we describe the data sets used for 
this analysis. We present our calculations of the quasar- 
galaxy correlation function in section 3. A discussion 
of the findings and analysis is presented in Section 4 & 
5 and in Section 6 we discuss our results and we give 
some concluding remarks. Throughout this paper we use 
cosmological parameters Q m = 0.3, ft a = 0.7 and Ho = 
70 km s _1 Mpc -1 . 

2. THE DATA 

2.1. The 2SLAQ Luminous Red Galaxy Sample 

Our LRG sample is taken from the 2SLAQ survey de- 
scribed in Cannon et al. 2006 (hereafter C06). The pub- 
lic data set (2SLAQ LRG Redshift Catalogue, version 5 
available at www.2slaq.info) contains nearly 15000 LRGs 
with magnitudes brighter than i = 19.8 along with their 
measured spectroscopic redshifts. We select only those 
regions from the northern stripe that overlap the 2QZ 
sample and do not use parts of the survey where the 
completeness is low (i.e., fields aOl and a02). We select 
galaxies that cover the redshift range from 0.2 to 0.8 with 
a median redshift of 0.52. Figure 1 shows a histogram of 



Fig. 2. — The figure shows the projected two point correlation 
function for 2QZ quasars with 2SLAQ LRGs plotted as diamonds. 
la 'Jackknife' errors are also shown. The best fit function of the 
form in Equation 4 is overplot (solid) for scales 0.7-27/i — ^^Mpc. 
Parameter values are ro = 5.3 ± 0.6/i~ x Mpc and 7 = 1.83 ± 0.42. 
The abscissa positions have been slightly offset for clarity. We also 
show our measured LRG-LRG autocorrelation function (squares) 
with the functional fit of r = T-l^o^h- 1 Mpc and 7 = 1.98±g;gg 
(dash-dot). 

galaxy redshifts. Our sample is selected from the 'top' 
priority sample described in C06 (Sample 8, selected with 
priority parameter=8). These are sources that fall in the 
most conservative r — i,g — r color cut of the 2SLAQ 
survey. We further selected only those sources with good 
spectral quality (Q > 3). The overall measured redshift 
completeness for the 2SLAQ Sample 8 (that is, the per- 
centage of the sample with good enough spectra taken 
to provide reliable redshifts) is 91.4% (C06, Ross et al. 
2007). The color and magnitude restrictions ensure that 
our full galaxy sample is dominated by the most massive 
and luminous intermediate redshift galaxies, these are 
expected to be galaxies with a passively-evolving stellar 
population (see Eisenstein et al. 2001). Our final galaxy 
sample contains 4985 LRGS. 

2.2. The 2dF Quasar Sample 

Although the SLAQ survey also includes quasars, we 
choose our quasar sample from the 2dF QSO Redshift 
Survey (2QZ, Croom et al. 2004). The primary reason 
for this is to avoid sampling problems resulting from sur- 
vey priorities, fiber placement, etc. within the single sur- 
vey. For the 2SLAQ sample, quasars are assigned lower 
priority in the placement of spectral fibers. Although 
there is overlap of the fields that would allow lower pri- 
ority sources to be picked up in subsequent observations, 
the subtle vagaries of target selection and fiber placement 
will lead to a bias in the measurement of the two point 
correlation function, particularly at smaller scales. In 
order to avoid having to make significant corrections for 
this 'crosstalk' between quasars and galaxies that may 
vary from field to field, we instead use the overlapping 
2QZ survey. This selection limits us to the study of more 
luminous quasars since the 2QZ targets are brighter than 
the 2SLAQ survey targets. Although there is no crosstalk 
between our two samples, there are still issues of fiber 
placement in both samples, thus we are limited to mea- 
suring correlation functions on intermediate scales be- 
tween about 0.7 - 30 Mpc. 

The 2QZ survey is a flux limited, optically selected 
quasar survey conducted within the 2df galaxy redshift 
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survey. The survey includes 23,338 quasars with bj mag- 
nitude range of 18.25 < bj < 20.85. Quasar candi- 
dates were selected with U-excess in the U — bj, bj — R 
color plane and confirmed with spectra of emission lines. 
We select QSOs that overlap the non-contiguous region 
of the 2SLAQ survey within the galaxy redshift range 
of 0.1 < z < 0.8. Figure 1 shows the redshift dis- 
tribution of our final quasar sample. Within this red- 
shift range, we have selected luminous quasars with 
-21 > M b j > -24 (Croom et al. 2004). From the pub- 
lic database (www.2dfquasar.org/Spec_Cat/), we select 
only those sources classified as 'QSO' that have quality 
flags of 11, indicating that the spectra are of high quality 
and redshift measurements are secure. Our final sample 
contains 420 quasars. 

3. MEASURING THE CROSS CORRELATION OF QUASARS 
AND LRGS 

The cross correlation function is used to measure the 
clustering of sources above a random distribution de- 
scribed by Poisson statistics. This correlation function 
is a measure of the excess probability of finding a source 
from one sample within a volume element, dV, at a dis- 
tance r from a source in a second sample. Here, we use 
a standard estimator (Landy & Szalay 1993) to mea- 
sure the quasar-LRG cross correlation function, £, where 
we measure the observed number of LRGs around each 
quasar as a function of distance and divide by the ex- 
pected number of galaxies from a random distribution. 



tl \ n Random QG(r) 
f (r) = * 



Tlgala 



QR(r) 



(1) 



Here QG is the number of quasar-galaxy pairs and QR is 
the number of quasar-random pairs at a given separation. 
n ra ndom/n galaxies is a normalization factor for the mean 
number densities in each catalog. 

We create random catalogs of LRGs that mimic the 
spatial and redshift distribution of our galaxy sample. 
Because of the non-contiguous nature of the 2SLAQ sur- 
vey in right ascension, in order to properly distribute the 
random sample galaxies, we use the same RA positions 
of the real galaxies. We then randomize the redshift and 
declination of each mock galaxy. For these catalogs, each 
source redshift (declination) is selected at random from 
the observed redshift (declination) probability distribu- 
tion. We use a density of random sources that is 50 times 
our galaxy density. 

We expect apparent distortions of the clustering pat- 
tern ('Fingers of God') due to peculiar motions of galax- 
ies along the line of sight (Peebles, 1980). To account for 
these distortions, we measure £ as a function of both r p , 
component perpendicular to the line of sight and ir, the 
component along the line of sight. The real space cor- 
relation £(r), without distortions, can then be recovered 
by integrating £(r p , n) along the line of sight to calculate 
co(r p ), the projected correlation function. 



uj{r p ) = 2 £(r p ,7r)d7r 
Jo 



(2) 



It has been shown by Davis and Peebles (1983) that the 
projected correlation function is related to the real-space 
correlation function through the equation, 
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Fig. 3. — In this figure, we compare our measured two point 
correlation function of 2QZ quasars with the full sample of 2SLAQ 
LRGs (diamonds) to those of other recent work. Mountrichas et 
al. (2008) find r = 8.68 ± 1.0 and 7 = 1.63 ± 0.15 over scales of 1 
to 25 h~ Mpc. Their selection criteria is somewhat different than 
ours limiting the redshift range to 0.35 < z < 0.75 and effectively 
eliminates more bluer galaxies. There result gives a much larger 
ro and and a more shallow slope. We also plot the correlation 
function from Padmanabhan ct al. (2008) between SDSS quasars 
and LRGs at 0.2 < z < 0.6. They measure a correlation function 
fit by a power law slope of 1.8 ± 0.1 and ro = 6 ± 0.5h -1 Mpc over 
scales of 0.1-20/i _1 Mpc, which compares well to what we find. 
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Lo{r p ) = 2 / rdr£{r)(r 2 - r z p ) 
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(3) 



In practice, we integrate not to infinity but to some 
large value, 7r ma ~ — 50ft. _1 Mpc, which includes almost all 
correlated pairs and peculiar velocities. We have varied 
the value of ir m ax using both 30/i _1 and 50/i _1 Mpc and 
find that our results for r and 7 are changed by < 10%. 
Our amplitudes (measured at 1.7/i _1 Mpc) change by as 
much as 20%, however this change has no effect on our 
conclusions since amplitude differences between our two 
subsamples (see section 5.1) are larger than these values. 

The real-space correlation, £(r), can then be approxi- 
mated by a power law. The projected correlation func- 
tion is 



Lu{r p ) 



r(i/2)r(( 7 -i)/2) 
r(7/2)] 



(4) 



where T is the Gamma function. However, real-space 
correlation function is not expected to be a true power 
law on large scales (ro > 20/i _1 Mpc). Below we present 
the results of our cross correlation along with power law 
fits to the function. 

4. QUASAR-LRG CORRELATION RESULTS FOR THE 
FULL SAMPLE 

Figure 2 shows the measured projected correlation 
function between our selected sample of 2QZ quasars and 
2SLAQ LRGs. Although we have used the largest sample 
of LRGs with spectrally measured redshifts, our ability 
to measure the correlation on small scales (< 0.5Mpc) is 
hampered by the small number of quasar-galaxy pairs in 
each redshift bin. Errors have been estimated by using 
the 'jackknife' method of Scranton et al. (2002), where 
we divide the sample into 12 approximately equal area 
regions and calculate the correlation function using com- 
binations of only 11 regions, leaving out the one region 
for each re-calculation. 
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Fig. 4. — The top plot in this figure shows a histogram of rcdshifts 
for the red (2569) and blue (2416) 2SLAQ LRG samples partitioned 
at g — r = 1.6. The distributions of the redshifts are similar with 
medians of 0.54 and 0.55, respectively. The bottom plot is a his- 
togram of r magnitudes for the same red and blue 2SLAQ galaxy 
samples. Again, the distributions are very similar. 



DRjjrp) 
f~[ DR(r p ) 



[uJi{r p ) - uj(r p )} 2 



(5) 



DRi is the number of data-random pairs for each sub- 
sample, DR is the number of pairs in the full sample, 
thus the ratio is a normalization that accounts for the 
fact than the number of random pairs may not be the 
same for each subfield (Myers et al. 2005). 

We have done a \ 2 n ^ of the power law model de- 
fined in Equation 4 to the correlation over the range 
of 0.7 — 27/i~ 1 Mpc and find a clustering radius, r = 
5.3 ± 0.6 and slope, 7 = 1.83 ± 0.42. In Figure 3, we 
compare our results to Mountrichas, et al. (2008) who 
have also recently measured the cross correlation of 2QZ 
quasars and 2SLAQ LRGs. They find r = 8.68 ± 1.0 
and 7 = 1.63 ±0.15 over scales of 1 to 25 h~ 1 Mpc. They 
have used a different selection criteria to ours limiting 
the redshift range to 0.35 < z < 0.75. This selection 
effectively eliminates more bluer galaxies than red from 
the LRG sample at both ends of the redshift range (see 
figure 3). Thus they have selected a more intrinsically 
bright sub-sample of galaxies, irrespective of color (see 
Figure 4). If we use a selection more similar to theirs, 
we find a somewhat larger r , but we do not match their 
value. 

We have also measured the autocorrelation function 
for the sub-sample of 2SLAQ galaxies used in our se- 
lected sample. The auto correlation function is also plot- 
ted in Figure 2. A fit to the measured function gives 
r = 7.lt£;^- 1 Mpc and 7 = 1.98±q ^ for scales of 
0.7 to 47/i _1 Mpc. This result closely matches that of 
Ross et al.(2007) who find r = 7.30 ± 0.34/i" 1 Mpc and 
7 = 1.83 ±0.05 over similar scales from 0.4 to 50/i _1 Mpc 
from a measure of the projected correlation function of 
the full 2SLAQ galaxy sample. We therefore conclude 
that galaxy distribution properties or our selected galaxy 
sample are representative of the full 2SLAQ LRG sample. 

We confirm that the quasar-LRG correlation has a 
somewhat lower amplitude than the LRG-LRG correla- 
tion. In particular at scales < 3/i _1 Mpc, the COrrela- 
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Fig. 5. — The plot in this figure shows histograms of r-band 
absolute magnitudes for the red and blue LRG sample partitioned 
at g — r = 1.6. Shown here is that the red sample of galaxies is, on 
average, more luminous than the blue sample. We interpret this to 
mean that the red sample is more massive than the blue sample. 

tion function seems to be more shallow indicating that 
a single sloped power law may not be the best fit. We 
conclude that our luminous quasars are are strongly cor- 
related with these LRG galaxies, although they seem to 
be less strongly correlated at scales < 3/i _1 Mpc. To try 
and make some sense of what this means for quasar trig- 
gering scenarios, we take a look at the makeup of our 
LRG sample in the sections below. 

5. UNDERSTANDING THE LRG SAMPLE 

Roseboom et al. (2006) have looked in detail at the 
galactic content of the 2SLAQ LRG sample. Based on 
the equivalent widths of H<5 and [O/I] , they identify four 
spectral types of galaxies represented. The majority of 
the galaxies in the sample (81%) are found to be passively 
evolving, based on the absence of absorption or emission 
lines. Indeed, an inspection of the original selection of the 
2SLAQ LRG targets from Eisenstein et al. (2001) (also 
see Roseboom et al. 2006), shows that the evolutionary 
track model for a passively evolving early-type galaxy 
runs nearly vertical in color-color space at g — r ~ 1.7 
through the r — i range where these galaxies have been 
selected. 13% of the galaxies show [O II] emission with- 
out H5 absorption, while 4% show H5 absorption without 
significant [O II] and 2% have both H5 absorption and 
[Oil] emission. The presence of emission and absorption 
lines are believed to indicate ongoing activity in these 
galaxies; US absorption in 6% of these galaxies indicates 
recent starformation which makes these galaxies appear 
bluer in the LRG sample. Roseboom et al. (2006) sug- 
gest that this starformation is unlikely to be continuous 
since evolutionary tracks for such a scenario barely skim 
the color cut imposed by the 2SLAQ survey. Thus, they 
conclude that the starformation in these galaxies is the 
result of 'instantaneous' recent bursts. 

The sources in the 2SLAQ LRG sample that are classi- 
fied as emission line galaxies (~ 13%) may be the result 
of either ongoing starformation or AGN activity. Rose- 
boom et al. 2006 were unable to distinguish between 
these because of the lack of spectral coverage of the di- 
agnostic lines. Therefore, observationally ~ 20% of the 
2SLAQ LRG sample are seen to have recent active star- 
formation or perhaps even ongoing AGN activity. 

Wake et al. (2006) demonstrate that at the bluer end 
of the 2SLAQ sample, magnitude errors push less intrin- 
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Fig. 6. — Shown is the projected two-point correlation function 
for 2QZ quasars and blue galaxy sample with 1 a errors. The best 
fit to Equation 4 gives a ro = 7.7±0.8 and 7 = 2.17^q'3q for scales 
l-10h _2 Mpc. This is a stronger correlation than that for the full 
LRG sample. 
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Fig. 7. — The projected two-point correlation function for 2QZ 
quasars and red LRGs is shown in this plot. The correlation func- 
tion is noisy despite the larger number of galaxies in the red sam- 
ple. The fit for scales of l-10/i _2 Mpc give a ro = 0.9 ± 0.1 and 
7 = 1.07 ± 0.21, however the correlation flattens on scales < 5. 

sically luminous galaxies into the selected LRG sample. 
Their analysis of the photometric errors shows this to be 
a 22% effect in the full 2SLAQ sample. They also show 
(in Figure 1 of their paper) that for the redshift range 
covered, evolutionary models that include starformation 
are much better fits to the colors of galaxies blueward of 
g — r = 1.6. Thus our LRG sample seems to be made up 
of 2 different galaxy populations that have different SEDs 
and have had very different evolutionary histories. One 
less luminous blue population with on going activity and 
a second red, more luminous population that is passively 
evolving. In the next section, we examine the correlation 
of 2QZ quasars with these two different galaxy popula- 
tions. 

5.1. Quasar-LRG Correlation Results for Red/Blue 
sub-Samples 

A division of our selected sample at g — r = 1.6 leaves 
us with similar numbers of galaxies in the red (2569) and 
blue (2416) subsamples. Figure 3a shows a histogram of 
the redshift ranges for these red and blue subsamples for 
the top priority 2SLAQ sample and demonstrates that 
the distributions are quite similar with median redshifts 
of 0.54 and 0.55, respectively, suggesting that differences 
in the quasar-galaxy correlation function for each of these 



samples is not due to redshift evolution. A histogram of 
the r magnitudes for each of the samples is shown in 
Figure 3b. Again, these distributions are similar. 

Although these 2 galaxy sub-samples have similar dis- 
tributions in redshift and magnitude, we have shown 
above that they likely have different spectral energy dis- 
tributions (SEDs) - about half of galaxies in our blue 
sample are starforming, and possibly disk galaxies. Thus 
the absolute r-band magnitude distributions are not the 
same for these populations. We estimate the absolute 
magnitude differences in the distributions by assuming 
an early type galaxy SED for our red population and 
a late type distribution for our blue galaxies, which re- 
flects the increased star formation in these sources. We 
determine r-band K-corrections for these populations as- 
suming an average redshift of 0.55. Figure 4 shows the 
difference in the derived absolute magnitudes for each 
sample. The red sample is shown to have a higher av- 
erage luminosity than the the blue sample, as shown by 
Wake et al. (2006). We assume that mass follows light 
in these galaxies and that this higher luminosity implies, 
on average, higher masses for our red galaxies. 

As we did for the full LRG sample, we now measure the 
projected correlation function, ui p , for our red and blue 
galaxy subsamples with the full 2QZ quasar sample. The 
correlations are plotted in Figures 5 and 6 with their one 
sigma Poisson errors. We use the Poisson errors here 
because we found that the jackknife errors seem to mis- 
represent the true error in our correlation function. The 
size of the subsamples are not large enough to partition 
the samples into enough bins for the jackknife estimate to 
be robust while keeping the number statistics reasonable 
for each partition. Other studies (e.g. Ross et al. 2007, 
da Angela 2008) have shown that at least for small scales 
(< 10ft- _1 Mpc), where bins are independent, the Poisson 
errors are similar to the jackknife and field-to- field errors. 

The fitted ro and 7 are reported in Table 1 along 
with the reduced x 2 of the fit. We find that the cor- 
relation function for the quasar-blue galaxy sample on 
scales 0.7 to 27/i" 2 Mpc has a best fit r = 7.3 ± 0.7 
and 7 = 2.11 j[ 13 with higher amplitude than the cross 
correlation with the entire galaxy sample although the 
reduced \ 2 value indicates that correlation function is 
poorly fit by a single power law over this range. On scales 
of ~ 1 — 10ft~ 2 Mpc however the function is well fit by 
a power law with r a = 7.7 ± 0.8/i~ 2 and 7 = 2.17±g;^, 
an r value larger than that measured for the full galaxy 
sample. 

In contrast, the red sample correlation is very noisy, 
despite the larger numbers of galaxies in the sample. The 
power law fit over scales 0.7 to 27ft~ 2 Mpc gives a r = 
4.9 ± 0.7 and 7 = 1.66 ± 0.10 , but once again, these 
are poor fits to a single power law model. The fit over a 
narrower range of scales, ~ 1 — 10ft.~ 2 Mpc, gives a nearly 
flat fit with a formal r = 0.9 ± 0.1 and 7 = 1.07 ± 0.21. 
The data are very noisy because the correlation signal is 
weak, particularly on small scales. This is reflected in 
the red sample errors. 

We also measure the auto-correlation of our red and 
blue LRG samples, shown in Figure 7. Results of the 
power law fits are listed in Table 1. We find that the 
redder, more massive LRGs, have a higher clustering 
amplitude than the bluer galaxies. That is, there is a 
higher probability of finding a red-red galaxy pair than 
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TABLE 1 

Table of fit data for correlation functions. 
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h^Mpc 



X 2 

reduced 
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Fig. 8. — The projected auto correlation function for blue and 
red galaxies are shown along with the fit over scales 0.7-27/i" 2 Mpc. 
The blue galaxy data points have been slightly offset for clarity. 
Fitted values for ro and 7 are listed in table one. The red galaxies 
are more strongly clustered particularly on small scales. 



a blue-blue pair, at least at small scales. Thus the red- 
der galaxies are more strongly clustered than the bluer 
galaxies, although the clustering radius (where the dis- 
tribution becomes random) is larger for red galaxies. 

6. DISCUSSION 

From the above results, we conclude that most of the 
signal in the measured projected two point correlation 
of bright 2QZ quasars cluster tothe full sample of LRGs 
on scales < 5 Mpc is the result of a stronger correlation 
of quasars with the bluer galaxies in the LRG sample. 
While there is a measured correlation with the red galax- 
ies, the signal is weak and flattens on scales < 5/i _1 Mpc 
(see Figure 5). 

Our result that 2QZ quasars seem to cluster more 
strongly with the bluer galaxies than red galaxies is an 
interesting one. The expected merger driven evolution- 
ary history of the reddest LRG galaxies suggests a sce- 
nario in which these galaxies occupy the highest mass 
halos where they are built up through major mergers of 
galaxies in the highest density regions, which trigger a 
luminous quasar phase, we might have expected to de- 
tect a strong correlation between the most luminous red 
galaxies and our bright quasars, indicating that they oc- 
cupy dark matter halos of similar mass. The lack of a 
strong correlation with bright quasars suggests that these 
quasars are not just found in the highest density regions. 
Our finding is qualitatively supported by Mountrichas et 
al. (2008), who find that luminous quasars may be less 
highly biased than faint quasars when compared with 
LRGs, that is, that the more luminous quasars arc not 
only hosted in the most overdense mass peaks. We ex- 
plore 2 explanations for our measured LRG subsamplc 
correlation functions. 

We have demonstrated that a large fraction of the blue 
galaxies in the LRG sample are generally lower lumi- 
nosity (~ L*) galaxies that have undergone some recent 
burst of starformation and a significant number of them 
may even host (low luminosity) AGN. Our strong corre- 
lation of quasars with these blue galaxies suggests that it 
is these galaxies that have space densities similar to our 
bright quasar sample. The fact that we exclusively mea- 
sure a strong correlation of luminous quasars with the 
population of galaxies that appear to have had recent 
activity (either an instantaneous burst of starformation 



or low luminosity AGN) supports models that predict 
the association of starbursts with AGN activity (e.g., 
Springcl et al. 2005). We suggest that environmentally 
dependent triggering mechanisms for the bursts of star- 
formation in these galaxies and the quasar activity may 
be the same. We are unable to say what that mechanism 
might be; that is whether minor mergers and harassment 
through tidal interactions might trigger activity or ma- 
jor merger events are contributing. However, if some of 
these lower luminosity galaxies do have disk components, 
they are not expected to be the resultant population of 
major mergers in hierarchical models. 

This suggestion is strengthened by the observation that 
red galaxies in our sample do not have space densities 
comparable to or bright quasar distribution. Our mea- 
surement of a turnover in the correlation function on 
scales < 5/i _1 Mpc between quasars and only the red, 
more centrally concentrated galaxies reflects the lack of 
quasar-red LRGs pairs above a random distribution on 
small scales. This can be explained if a significant frac- 
tion of 2QZ quasars are distributed in regions outside 
the more centrally clustered red galaxies. This leads us 
to conclude that the processes responsible for growing 
LRGs (presumably major mergers) are not those that 
responsible for triggering most bright quasars at this red- 
shift. 

Recent work on the growth of galaxies in hierarchical 
simulations by Guo and White (2008) lends support to 
our conclusions that major mergers may not be the en- 
tire story when it comes to triggering luminous quasar. 
They look at the growth of galaxies through major merg- 
ers, minor mergers and starformation for galaxies in a 
series of binned mass ranges over cosmic time. They find 
that the dominant mode of galaxy assembly is a strong 
function of the stellar mass of the galaxies. For galaxies 
with lower masses (M < 8 x 1O 1O M ), starformation 
and minor mergers are the dominant modes for galaxy 
growth. Whereas for higher mass galaxies major mergers 
are much more important for growing the galaxies at all 
rcdshifts. This work mirrors what has been found in the 
2SLAQ sample, that the lower mass (blue, L») galaxies 
are more closely associated with starformation activity, 
here in the form of bursts, than the more massive (red) 
galaxies. The strong correlation of quasars with these 
blue galaxies, suggests that minor mergers should be im- 
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portant triggers of activity, even for these bright quasars. 

However, Wake et al. (2008) conclude that there is ma- 
jor merging going on within the LRGs population over 
rcdshifts of 0.5 to 0.19. By comparing the autocorre- 
lation function at high and low-z with passive evolution 
and HOD models, they find that the merger rate over this 
redshift range is 2.4% per Gyr. The 2SLAQ LRG sam- 
ple they use is most similar to our red galaxy sub-sample. 
Thus for an assumed scenario where major mergers result 
in quasars, we would expect a strong correlation of red 
LRGs and quasars. So why is this not what we measure? 
It is likely that our selection of bright 2QZ quasars does 
not enable us to measure a correlation because quasars 
resulting from these mergers arc faint due to dust obscu- 
ration. This is consistent with work suggesting that faint 
quasars are more highly biased than bright quasars. 

We also point out the recent work of Padmanabhan et 
al. (2008) on correlations of SDSS quasars and LRGs 
at 0.2 < z < 0.6. Although they have a larger number 
of LRGs (chosen similarly to our sample) and quasars 
in their sample, they do not have the directly measured 
rcdshifts for each galaxy, as we have. They measure a 
correlation function fit by a power law slope of 1.8 ± 0.1 
and ro = 6±0.5/i~ 1 Mpc over scales of 0.1-20/i _1 Mpc (see 
Figure 3). With their very large numbers of quasars and 
galaxies, they are able to measure the correlation func- 
tion to smaller scales than is possible for our sample. 
They model their correlation function with halo occupa- 
tion distribution models where quasar host galaxies pop- 
ulate their dark matter halos preferentially at the centers 
and with varying fractions of quasar hosts distributed as 
'satellite' galaxies within halos. They find that in order 
to match the correlation function, a significant fraction 
(> 25%) of quasars must be hosted outside the central 
halo region in satellite galaxies. 

Coil et al. (2007) have looked at the correlation func- 
tion of very luminous (-25.2 < M bj < -22.1) SDSS 
quasars with galaxies in the DEEP2 field at about red- 
shift = 1.0. They find that quasars are found in re- 
gions with mean overdensities more similar to blue galax- 
ies than red galaxies. We note however that their blue 
galaxy sample is, on average, much more blue than our 
blue sample. Their selection of blue galaxies are those 
that have U — B < 1 and constitute a more clearly disk 
dominated population. Because targets for the 2SLAQ 
survey were chosen to mostly exclude such galaxies, we 
arc unable to compare directly to a similar sample at 
lower redshift here. However, although a smaller fraction 
of our galaxies might be disk dominated, the similarity 
is in what makes these galaxies blue and that is the in- 
crease in starformation over the reddest LRGs. Like the 
higher redshift sample, we too find that quasars at z=0.5 
are closely associated and cluster like galaxies with en- 
hanced starformation. 

7. SUMMARY AND CONCLUSIONS 

We have measured the projected two-point cluster- 
ing of quasars from the 2QZ survey and 2SLAQ LRG 
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galaxies. We find that ro = 5.3 ± 0.6ft- 2 Mpc and 
7 = 1.83 ± 0.42 on scales 0.7 to 27 /i^Mpc. A liter- 
ature search of the galaxy populations that make up this 
full LRG sample, revealed that our galaxy sample con- 
tains two galaxy populations that might have very dif- 
ferent evolutionary histories. Thus we sought insight to 
our goal of exploring quasar triggering by looking at the 
clustering signal with each of these galaxy populations. 

We find that the clustering signal measured with the 
full galaxy sample can be attributed primarily to the 
strong clustering of our quasar sample with the bluer 
LRGs. Analysis of the spectra of these blue galaxies 
taken from the literature shows that they are primarily 
galaxies that have most likely undergone some recent in- 
stantaneous burst of starformation or may even harbor 
low luminosity AGN. 

We further find that the correlation signal with the 
sample of the most luminous and red LRGs is weak, 
noisy, and seems to flatten on scales < 5 Mpc. This 
strong correlation with blue galaxies and flattening with 
red, more centrally clustered galaxies can be explained 
by either a significant fraction of quasars being hosted in 
galaxies outside the cores of clusters or/and that bright 
quasars are less highly biased than faint quasars so that 
our selection criteria excludes most of the faint quasars 
that are correlated with our red sample of LRGs. 

We speculate on what this correlation means for the 
triggering of luminous quasars. The strong correlation 
between bright quasars and the bluer galaxies in our sam- 
ple LRG suggests that the triggering of these quasars is 
related to the starformation in these galaxies at z=0.5. 
Other studies show that this activity is more likely to 
be the result of minor mergers or tidal effects between 
galaxies than major merger events for these lower mass 
galaxies. Thus minor mergers and/or tidal effects may 
be important triggers of even bright quasars. 

If the spatial distribution of our quasars compared to 
LRGs is responsible for our observed stronger correla- 
tion with bluer galaxies, then this may support the so 
called "cosmic downsizing" scenario (Cowie 1996, Broad- 
hurst et al. 1988), where at higher redshift, quasar ac- 
tivity was more prevalent in high mass galaxies. By red- 
shift 0.5, this activity has ended, perhaps leaving behind 
a passively evolving LRG. The bulk of activity at this 
lower redshift is now in less luminous, lower mass galax- 
ies. However, if the redder LRGs that are still merging 
lead to quasar ignition, then there may still be some of 
quasar activity in high mass galaxies at z=0.5, but it is 
obscured by dust. 
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